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ABSTRACT 

Type Iln supernovae (SNe Iln) are a rare (< 10%) subclass of core-collapse SNe that exhibit 
relatively narrow emission lines from a dense, pre-existing circumstellar medium (CSM). In 2009, 
a warm Spitzer Space Telescope survey observed 30 SNe Iln discovered in 2003-2008 and detected 10 
SNe at distances out to 175 Mpc with unreported late-time infrared emission, in some cases more 
than 5 years post-discovery. For this single epoch of data, the warm-dust parameters suggest the 
presence of a radiative heating source consisting of optical and X-ray emission continuously generated 
by ongoing CSM interaction. Here we present multi-wavelength follow-up observations of this sample 
of 10 SNe Iln and the well-studied Type Iln SN 2010J1. A recent epoch of Spitzer observations 
reveals ongoing mid-infrared emission from nine of the SNe in this sample. We also detect three of 
the SNe in archival WISE data, in addition to SNe 1987A, 2004dj, and 2008iy. For at least five of 
the SNe in the sample, optical and/or X-ray emission confirms the presence of radiative emission 
from ongoing CSM interaction. The two Spitzer nondetections are consistent with the forward shock 
overrunning and destroying the dust shell, a result that places upper limits on the dust-shell size. 
The optical and infrared observations confirm the radiative heating model and constrain a number 
of model parameters, including progenitor mass-loss characteristics. All of the SNe in this sample 
experienced an outburst on the order of tens to hundreds of years prior to the SN explosion followed 
by periods of less intense mass loss. Although all evidence points to massive progenitors, the variation 
in the data highlights the diversity in SN Iln progenitor evolution. While these observations do not 
identify a particular progenitor system, they demonstrate that future, coordinated, multi-wavelength 
campaigns can constrain theoretical mass-loss models. 
Subject headings: circumstellar matter — supernovae: general — supernovae: individual (SN 2010J1, 

SN 2008gm, SN 2008en, SN 2008J, SN 2006jd, SN 2005ip) — dust, extinction - 

infrared: stars 



1. INTRODUCTION 

Type Iln supe rnovae (SNe Iln; iSchlegell I1990L see 
iFilippenkol 119971 for a review) have gained consider- 
able attention over the past decade given their associa- 
tion with explosions r anging from underluminous "super- 
nova impostors" fe.g.JSmith et al.ll2011l ; lKochanekll2011l ; 



yr _1 ) to episodic dense winds observed in some 
massive stars (~20-40 M or greater ), including a 
small fraction of red supe rgiants (RSGs; iFransson et al.l 
120021 ISmith et al.l 12009a ) and luminous blue variables 



(LBVs; e.g.. ISmith et al. 



120111 and references therein) . 



iVan Dvk fe Matheson I2012D to supe r luminous Type I I 
SNe (SLSNe: e.g.. ISmith et al.l [20071 IQfek et al.l [2001 . 
Representing fewer than 10% of all core -collapse event s 
in the nearby Universe (d < 60 Mpc; iLi et al.l 1201 ID . 
SNe Iln are characterized by relatively narrow lines (~ 
100-1000 km s -1 : lSchlegeilll990f ). This defining charac- 
teristic is not associated with the SN explosion itself, but 
rather with a dense circumstellar medium (CSM) pro- 
duced by pre-SN mass loss (though in some cases much 
of the line width can be attributed to interaction with 
the SN ejecta). 

Identification of a single progenitor class remains am- 
biguous. A popular viewpoint connects some high 
mass-loss properties (e.g., M s=s 10 _1 to 10 -3 M Q 
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IGal-Yam et all (120071) and IGal-Yam fc Leonard! (|2CJ09h 
also show direct observational evidence for a LBV pro- 
genitor in the case of SN 2005gl. 

The high mass loss, however, is not typical of most 
LBVs, which achieve rat es of only M < 10~ 4 Mp, yr" 1 i n 
their S Doradus state (Hump hreys fc Davidson! 119941 ). 
Furthermore, stellar evolution models are inconsistent 
with LBV progenitors. Classical stellar evolution the- 
ory doesn't even allow for terminal explosions in the 
LBV state, r equiring them to first evolve into a Wolf- 
Rayet phase (| Abbott fc Contil H9871 ) . The models also 
have difficulty achieving such substantial mass loss. At 
best, the mass-loss calculations rely on observed rates o r 
scaled line-driven wind models (jSmith fc Owockill2006|) . 
Alth ough some recent scenario s present promising leads 
(e.g.. lQuataert fc S hiodc 20l3), these models are incom- 
plete and cannot explain why only a fraction of massive 
stars undergo such extreme mass loss. 

While the exact mass-loss mechanism remains un- 
known, the pre-SN mass-loss history offers an important 
constraint on different progenitor models since varying 
wind speeds, densities, and asymmetries result in differ- 
ent observational behaviors. For example, a dense and 



Table 1 

Spitzer Observations 1 



SN 


JD 


Epoch 


a (h:m:s) 


8 (»:':") 


Distance 2 


'hit 


3.6 /im 3 


4.5 /im 3 




-2,450,000 


(days) 


(J2000) 


(J2000) 


(Mpc) 


(s) 


(10 17 erg s" 


1 cm- 2 A- 1 ) 


2002bu 


5734 


3372 


12:17:37.18 


+45:38:47.4 


6.9 


600 


<0.11 


<0.09 


20031o 


5865 


2860 


03:37:05.12 


-05:02:17.3 


57 


600 


<0.11 


<0.09 


2005cp 


5797 


2254 


23:59:30.88 


+18:12:09.6 


91 


300 


0.23(0.07) 


0.15(0.05) 


2005gn 


5875 


2219 


05:48:49.07 


-24:22:45.5 


165 


600 


<0.11 


<0.09 


2005ip 


5737 


2057 


09:32:06.42 


+08:26:44.4 


31 


300 


7.27(0.27) 


7.64(0.43) 


2006jd 


5735 


1714 


08:02:07.43 


+00:48:31.5 


77 


300 


2.68(0.16) 


2.66(0.24) 


2006qq 


5875 


1805 


05:19:50.43 


-20:58:06.4 


119 


300 


0.44(0.10) 


0.28(0.05) 


2007rt 


5733 


1304 


11:02:34.29 


+50:34:58.5 


93 


300 


1.63(0.14) 


1.52(0.12) 


2008J 


5842 


1361 


02:34:24.20 


-10:50:38.5 


66 


300 


<0.11 


<0.09 


2008cg 


5818 


1226 


15:54:15.15 


+10:58:25.0 


151 


600 


0.11(0.04) 


0.10(0.03) 


2008en 


5813 


1131 


00:55:13.56 


+35:26:26.2 


151 


600 


0.21(0.05) 


0.19(0.04) 


2008gm 


5797 


1035 


23:14:12.39 


-02:46:52.4 


19 


300 


0.14 (0.05) 


0.07 (0.03) 


2010J1 


5733 


230 


09:42:53.33 


+09:29:41.8 


50 


600 


13.08(0.40) 


8.59(0.40) 



1 Upper limits for nondetections were derived by the point-source sensitivity in Table 2.10 of the IRAC 
Instrument Handbook, version 2. 

2 Aside from SN 2002bu, all distances are derived from the host-galaxy redshift assuming Hq = 72 km 
s _1 Mpc -1 . 

3 lcr uncertainties are given in parentheses. 



optically thick CSM wind will increase the r adius of shock 
breakout relative to the stellar s urface (e.g., Moriva et al.l 
[20111: iChevalier fc Irwinl |201ll) . lOfek et al.l (|2012[ ) show 



the relationship bet ween mass loss and the time scalc for 
the shock breakout. IMoriva fc Tominagal (|2012t) further 
use nonsteady mass loss to explain the spectral diversity 
of Type II SLSNe (i.e., those with H in their spectra). 
While several core-collapse SNe (CCSNe) not classified as 
SNe Iln h ave shown evidence for late-time CSM interac - 
tion (e.g.. lAndrews et al.ll20To1 : iMilisavlievic et al.ll2012f ). 
the dense CSM associated with SNe Iln makes them a 
particularly good sample for late-time studies (> 100 d) 
to trace the complete mass-loss history of the progenitor. 
Estimated to have a rate of no mor e than 10 yr -1 out 
to 150 Mpc (jDahlen fc Frans son 1999J) newly discovered 
SNe Iln are not easy targets. Instead, iFox et all (J2011I ) 
executed a warm Spitzer Space Telescope (Spitzer) sur- 
vey (P60122) that included 30 SNe Iln previously discov- 
ered in 2003-2 008. The survey , along with observations 
of SN 2005ip (IFox et al.l [20091 I27M . revealed 10 SNe 
with previously unreported late-time infrared (IR) ex- 
cesses, in some cases more than 5 yr post-discovery. The 
data are most consistent with pre-existing dust shells 
with radii >0.01 pc, formed by a progenitor outburst 
>10-100 yr prior to the SN. The outbursts had mass- 
loss rates of M « 10 _1 to 10~ 3 M Q yr -1 and large to- 
tal masses M to t ~ 1-10 M Q . The dust is most likely 
continuously heated by visible and X-ray radiation gen- 
erated by ongoing CSM interaction, although direct evi- 
dence for such shock interaction in th is sample exists only 
for SNe 2 05ip 2 006.id, and 2007r t (ISmith et al.ll2009bt 
Fox et all 120111: iStritzinger et al.l 120121: iTrundle et al.l 



2009^ 



X-ray and optical observations of the shock interaction 
not only validate this model but also probe the inner 

< 0.01 pc of the CSM, which corresponds to mass loss 

< 100 yr prior to the SN. This paper presents a multi- 
wavelength follow-up study of the 10 SNe observed by the 
initial Spitzer survey (as well as the well-studied Type 
Iln SN 2010J1). Data include an additional epoch of warm 
Spitzer photometry, photometry from the Wide-field In- 
frared Survey Explorer (WISE), and ground-based opti- 



cal photometry and spectroscopy. The spectral energy 
distribution (SED) provides a more complete picture of 
the evolution of SNe Iln as a forward shock continues to 
pass through a dense CSM. Section [2] lists the details of 
our observations. In $3]we analyze the data with respect 
to the interaction model. Implications on the progenitor 
evolution models are discussed in g] Section [S] summa- 
rizes the results. 

2. OBSERVATIONS 
2.1. Warm Spitzer /IRAC Photometry 

Warm Spitzer /IRAC (|Fazio et al.l 120041 ) obtained a 
second epoch of data for the 10 SNe Iln disco vered to 
have late-time IR emission bv IFox et al.l (|2011l ) and the 
well-studied Type Iln SNe 2002bu, 20031o, and 2010J1 
(Program PID 80023). Table [1] lists the observational 
details. The Spitzer Heritage Archive] provided access 
to the Post Basic Calibrated Data (pbed), which are al- 
ready fully coadded and calibrated. The background flux 
in most of the SN host galaxies is bright and exhibits 
rapid spatial variations. Although template subtraction 
is a commonly used technique to minimize photometric 
confusion from the underlying galaxy, no pre-SN observa- 
tions exist. Instead, photometry was performed using the 
DADPHDT point-spread function (PSF) photometry pack- 
age in IRAF (see Figure [Q) 

2.2. WISE Photometry 

WISE mappe d the entire sky du ring 2010 at 3.4, 4.6, 
12, and 22 ^m (|Wright et al.ll2010n . WISE therefore in- 
cluded the same SNe that were observed by Spitzer. The 
NASA/IPAC Infrared Science Archived provided fully 
calibrated images of each field. Like the Spitzer data 
described above, photometry was performed using the 



SHA 



be 



accessed 



from 



http://sha.ipac.caltcch.edu/applications/Spitzcr/SHA/ 

D IRAF: the Image Reduction and Analysis Facility is distributed 
by the National Optical Astronomy Observatory, which is oper- 
ated by the Association of Universities for Research in Astronomy 
(AURA) under cooperative agreement with the National Science 
Foundation (NSF). 

7 IRS A can be accessed from [http://irsa.ipac. caltech.edu| 
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Table 2 

WISE Photometry 1 



SN 



JD 
-2,450,000 



Epoch 

(days) 



Wl 2 
(10 17 erg s^- 1 cm 



-1 



W2 2 

2 A- 



l ) 



1987A 


5311 


8461 


3.12(0.03) 


2.98(0.02) 


2002bu 


5232 


2871 


<0.27 


<0.23 


20031o 


5345 


2341 


<0.27 


<0.23 


2004dj 


5285 


2068 


0.34(0.06) 


0.51(0.03) 


2005cp 


5375 


1833 


<0.27 


<0.23 


2005gn 


5267 


1612 


<0.27 


<0.23 


2005kd 


5258 


1572 


<0.27 


<0.23 


2006gy 


5238 


1242 


<0.27 


<0.23 


2006jd 


5303 


1283 


2.79(0.04) 


2.43(0.02) 


2006qq 


5244 


1175 


<0.27 


<0.23 


2007rt 


5325 


897 


5.98(0.25) 


2.66(0.07) 


2008J 


5400 


920 


<0.27 


<0.23 


2008en 


5396 


715 


0.37(0.03) 


0.23(0.02) 


2008gm 


5356 


595 


<0.27 


<0.23 


2008iy 


5244 


507 


2.42(0.04) 


1.42(0.03) 



1 Upper limits for nondetections (~0.11 mjy at 4.6 )im) wer e 
set by the WISE point-source sensitivity (Wright ct al .120101) . 

2 la uncertainties are given in parentheses. 



DADPHDT PSF photometry package in IRAF. Pixel fluxes 
were converted from MJy sr" 1 to mjy according to the 
Explanatory Supplement to the WISE All-Sky Data Re- 
lease Products, which discusses the pixel size, zero points, 
and aperture correction in detail. Poor resolution and 
sensitivity resulted in nondetections for all SNc in the 12 
and 22 /im bands, so we do not include that photometry 
here. Table [5] lists the fluxes and upper limits. 

2.3. Optical Photometry and Spectroscopy 

Tables [3] and 0] and Figures [5] and |3] summarize the 
optical photometry and spectra. Data were obtained 
with the d ual-arm Low Re solution Imaging Spectrome- 
ter (LRIS; lOke et al|[l995l ) mounted on the 10-m Keck I 
telescope and the DEep Imaging M ulti-Object Spectro- 
graph fDEIMOS: lFaber et alil2003D mounted on the 10- 
m Keck II telescope. Keck/LRIS spectra were obtained 
using the 600/4000 or 400/3400 grisms on the blue side 
and the 400/8500 grating on the red side, along with 
a 1" wide slit. This resulted in a wavelength cov- 
erage of 3200-9200 A and a typical resolution of 5- 
7 A. Kcck/DEIMOS spectra were obtained using the 
1200/7500 grating, along with a 1" wide slit. This re- 
sulted in a wavelength coverage of 4750-7400 A and a 
typical resolution of 3 A. Most observations had the 
slit aligned along the parallactic angl e to minimize dif- 
ferential light losses (|Filippenkolll982t) ; moreover, LRIS 
is equipped with an atmospheric dispersion corrector. 

The photometric images were reduced using standard 
CCD processing techniques in IDL and Python, utiliz- 
ing online astrometry programs SExtractor and SWarpQ 
(Perley 2012, private communication). Photometry was 
performed using the DAOPHOT PSF photometry package 
in IRAF. Calibration was performed using held stars with 
reported fluxes in the Slo an Digital Sky S urvey (SDSS) 
Data Release 9 Catalogue (|Ahn et al.ll2012[ ). The spectra 
were r educed using standard techniques fe.g.. lFolev et al.l 
120031 : iSilvermanet alll2012t ). Routine CCD processing 
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Figure 1. Photometry of the SNe detected by Spitzer /IRAC 
bands 1 (3.6 /im) and 2 (4.5 /im). Overplotted are the resulting 
best fits of Equation [T] 



and spectrum extraction were completed with IRAF, and 
the da t a wer e extracted with the optimal algorithm of 
iHornd (J1986D . We obtained the wavelength scale from 



low-order polynomial fits to calibration-lamp spectra. 
Small wavelength shifts were then applied to the data 
after cross-correlating a template sky to an extracted 
night-sky spectrum. Using our own IDL routines, we fit 
a spectrophotometric standard-star spectrum to the data 
in order to flux c alibrate the SN and to remove telluric 
absor ption lines (jWade fc Hornel U988; Mathcso n et al.l 
[20(51 . 

2.4. Dust Temperature and Mass 

Assuming only thermal emission, the mid-IR photome- 
try provides a strong constraint on the dust temperature, 
mass, and thus IR luminosity. Assuming optically thin 
dust with mass M d and particle radius a, at a distance 
d from the observer, thermally emitting at a single equi- 
librium temperat ure T d , the flux can be written as (e.g., 
IHildebrandl[l98l 



F v = 



M d B u (T d )K u (a) 
d 2 



(1) 



8 S Extractor and SWarp can 
|http://www.astromatic.net/software| 



be accessed from 



where B v (T d ) is the Planck blackbody function and K u (a) 
is the dust mass absorption coefficient. 

For simple dust populations of a single size composed 
entir ely of either silicat e or graphite, the IDL MPFIT func- 
tion (jMarkwardtl 120091 ) finds the best fit (see Figure [J) 
of Equation Q] by varying M d and T<j to minimize the 
value of x 2 . The absorption coefficients, k, arc given 
in Figure 4 of iFox et ail (|2010t ). With only two mid-IR 
fluxes, our fits are limited to a single component. Table 
[5] lists the best-fit parameters for graphite grains of size 
a = 0.1 /j,m (found to be the best-fit value for these SNe 
bv lFox et al.ll201ll ). Figure [2] plots the corresponding IR 
luminosity evolution for each SN, as well as the luminosi- 
ties derived from i?-band and X-ray observations. 

3. ANALYSIS 

IFox et all (|2011f ) show that the dust properties for 
most of the 10 SNc Iln in their sample arc consistent 
with a pre-existing dust shell formed in a progenitor 
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Figure 2. The multi- wavelength luminosity evolution of the SNe detected by Spitzer. For reference, the near-IR evolution of two of the 
most well-studied SNe Iln (1995N and 2005ip) are overplotted (dashed and dotted lines). The mid-IR luminosities originate from warm 
dust, while the optical and X-ray luminosities likely originate from CSM interaction. While there is a certain degree of scatter in this plot, 
there is also a noticeable consistency. The mid-IR luminosities typically fall in the range 10 8 < £ m id-lR < 10 9 Lq, while the optical and 
X-ray luminosities arc slightly lower. Described further in the text, these values are consistent with a warm dust shell continuously heated 
by optical and X-ray emission generated by ongoing CSM interaction. 

Table 3 

Keck/LRIS Optical Photometry 1 



SN 


JD 
-2,450,000 


Epoch 

(days) 


B 

(mag) 


V 

(mag) 


R 

(mag) 


(mag) 


Luminosity 
(log (L opt /L )) 


2005cp 


5775 
5921 
6217 


2233 
2379 
2675 


>22.7 
>21.7 
>22.9 


>22.7 
>22.4 


>21.7 
>22.4 
>22.4 


>21.6 
>21.6 


<6.7 
<6.4 
<6.4 


2006jd 


5629 
6046 


1609 
2026 


22.5 (0.03) 
>23.2 


21.8 (0.03) 
>22.0 


20.4 (0.01) 
>21.8 


21.3 (0.02) 
>21.0 


7.57 
<6.7 


2007rt 


5629 


1223 


>23.3 


>22.9 


>22.7 


>22.3 


<6.3 


2008cg 


5774 
6094 


1183 
1503 


>23.2 
>22.5 


>23.1 
>22.1 


>22.5 
>21.7 


>22.1 
>21.2 


<6.4 
<6.7 


2008cn 


5775 


1094 


22.1 (0.04) 




21.6 (0.03) 


20.9 (0.03) 


7.63 


2008gm 


5774 


1013 


21.8 (0.04) 


21.8 (0.03) 


21.6 (0.04) 


20.7 (0.05) 


6.65 



Upper limits for nondetections reported at the 3a limit. 



wind/eruption and continuously heated by visible, ul- 
traviolet, and/or X-ray emission generated by ongoing 
CSM interaction in the forward shock. This scenario 
yields two testable predictions: (1) Assuming a single- 
shell model (Ar/r = 1/10), the grains will be sputtered 
and destroyed when the forward shock reaches the pre- 
existing dust shell, resulting in a decline in the IR light 
curve. If the CSM density drops sufficiently, the degree 
of shock interaction and, thereby, the optical light curve 
will also decrease. A declining IR light curve constrains 
the dust-shell radius and thus the heating model. (More 
complex CSM g eometries, such as mu ltiple dust shells 
(e.g., SN 2009ip; lMa"uerhan et all!2012[) . may complicate 



such an analysis.) (2) Assuming an optically thin dust 
shell, the observed dust temperature, T&, and shell ra- 
dius, r<j, require an optical, ultraviolet, and/or X-ray flux 
given by 



G4 



J opt/X 



par^aT^ 



SN 



J B v (T d )K{u)dv 

/S I/ (T SN )<5 a b s (^)d^ 



(2) 



for a dust bulk (volume) density p and an effective SN 
blackbody temperature Tsn, where B u is the Planck 
blackbody function, <5 a b s is the dust absorption effi- 
ciency, and k{v) is the dust mass absorption coefficient. 
Here we test these predictions. 
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Table 4 

Summary of Optical Spectra 



SN 


JD 
-2,450,000 


Epoch 
(days) 


Instrument 


Blue 
Res. (A) 


Red 
Res. (A) 


Int 
to 


2005ip 


5921 
6246 


2242 
2567 


Kcck/LRIS 
Kcck/DEIMOS 


6.5 


7 
3 


900 
2400 


2005cp 


6246 


2704 


Keck/DEIMOS 




3 


2400 


2006jd 


5896 
6246 


1876 
2226 


Kcck/LRIS 
Keck/DEIMOS 


6.5 


5.8 
3 


600 
2400 


2008en 


5866 
6193 


1185 
1512 


Kcck/LRIS 
Keck/DEIMOS 


3.6 


6.2 
3 


900 
2400 


2008gm 


5866 
6193 


1105 
1432 


Kcck/LRIS 
Keck/DEIMOS 


3.6 


6.2 
3 


900 
2400 



Late-Time Keck Spectra 




05ip, d 2567 
05cp, Id 2704 
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Figure 3. Optical spectra of SNe 2005cp, 2005ip, 2006jd, 2008en, and 2008gm taken with Keck/LRIS and DEIMOS. All spectra were 
obtained at > 3 yr post-explosion. While the spectra lack any significant continuum, they all show a significant Hq line, while SNe 2005ip 
also exhibits high-ionization Fe lines. 



Table 5 

IR Fitting Parameters (a = 0.1 /im) 



SN 


Telescope 


T d (K) 


M d (M Q ) 


L d (L ) 


2004dj 


WISE 


486 


2.39X10" 5 


9.19x10" 


2005cp 


Spitzer 


731 


5.17X10" 4 


2.20X10 7 


2005ip 


Spitzer 


511 


1.93X10" 2 


9.98X10 7 


2005kd 


WISE 


795 


4.23xl0~ 3 


2.87x10 s 


2006gy 


WISE 


1287 


5.57X10" 3 


6.52X10 9 


2006jd 


Spitzer 


531 


3.49X10" 2 


2.26x10 s 


2006jd 


WISE 


574 


2.08X10- 2 


2.11x10 s 


2006qq 


Spitzer 


771 


1.32xl0~ 3 


7.47X10 7 


2007rt 


Spitzer 


555 


2.25xl0~ 2 


1.87x10 s 


2007rt 


WISE 


834 


5.80xl0~ 3 


5.19x10 s 


2008cg 


Spitzer 


575 


3.14xl0~ 3 


3.24xl0 7 


2008cg 


WISE 


807 


2.74X10" 3 


2.02 x10 s 


2008cn 


Spitzer 


563 


6.77X10" 3 


6.17X10 7 


2008en 


WISE 


706 


2.63xl0~ 3 


8.89X10 7 


2008gm 


Spitzer 


1074 


1.61X10" 5 


6.41 xlO 6 


2008iy 


WISE 


730 


2.01X10" 2 


8.29x10 s 


2010J1 


Spitzer 


735 


8.47xl0~ 3 


3.62x10 s 



3.1. Dust-Shell Destruction 

The forward shock will eventually destroy the dust via 
sputtering, resulting in a declining mid-IR component. 
For a spherically symmetric dust shell, the shock will 
reach the dust shell at a time t > r^b/vs, where v s is the 
forward-shock velocity and r\^ is the blackbody radius, 
which is derived from the IR luminosity in Table [S] assum- 
ing an optically thick shell. This equation places a lower 
limit on the timescale over which dust destruction can 
occur. For example, the forward shock may decelerate 
with time. Also, an optically thin dust shell composed 
of smaller grains emitting with modified blackbodies re- 
quires a larger shell radius to generate the observed flux. 
Both of these effects increase the time over which it takes 
the forward shock to reach the dust shell. 

For each SN detected by Spitzer, Figure H] plots the 
blackbody radius and dust temperature. The diagonal 
arrows signify that these values are upper limits given the 
limitations of the single-component fit in Figure [T] While 
the numbers vary slightly, a blackbody radius r\,b ^S 0.1 
ly and a shock velocity v s w 3000 km s _1 will result in 
a mid-IR plateau that continues for roughly t < 10.5 yr. 
Once the shock reaches the shell, the luminosity decline 
will occur on timescales consistent with the dust-shell ra- 
dius (tdeciine = 2rbb/c ~ 0.2 yr), as the dust destruction 
on the near side of the shell will be observed before dust 
destruction on the far side. 

We search for evidence of dust destruction between the 
first and second epochs of Spitzer observations. Figure 
[2] shows that SNc 2005gn and 2008J are not detected 
on the second epoch of observations, which occurred at 
6.1 and 3.7 yr post-explosion, respectively. Assuming a 
3000 km s _1 shock, these results place upper limits on 
the dust-shell radii of r& <0.06 and 0.036 ly. A slower 
shock would result in a proportionally smaller shell ra- 
dius. While still detected, SN 2005cp also shows evidence 
of fading nearly 6.2 yr post-explosion. The horizontal 
arrows in Figure 0] illustrate both the lower and upper 
limits on the dust-shell radius derived for SNe 2005gn 
and 2008J. 
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Figure 4. Plot of the observed dust temperature, T^, and black- 
body radius, r^b- The diagonal arrows signify that these values 
reflect upper limits given the limitations of the single-component 
fit in Figure [T] The horizontal arrows for SNe 2005gn and 2008gn 
signify both the lower and upper limits on the dust-shell radius 
derived from the data. Overplotted are functions of the theoreti- 
cal dust temperature as a function of radius for contours of constant 
optical/X-ray luminosities. The results highlight that, in almost 
every case, the warm dust can be powered by optical and/or X- 



ray luminosities of 10 < L t 



pt/X 



< 10 Lq. Accounting for ab- 



sorption, these luminosities arc consiste nt with the observa tions 
shown in Figure [2] and those reported by [Smith et al. (2009^) and 
Chandra et al. (2012a|). 



3.2. Shock Interaction and CSM Characteristics 

For the scenario in which the dust shell is heated by 
continuous radiation generated by ongoing CSM inter- 
action, an optical and/or X-ray counterpart should ex- 
ist. Equation [2] calculates this luminosity as a func- 
tion of the observed dust temperature, Td, and shell 
radius, r<j. Along with the observed dust temperatures 
and radii, Figure H] plots the theoretical dust tempera- 
ture as a function of shell radius for contours of constant 
luminosity. In most cases, the observed dust parame- 
ters are consistent with optical and/or X-ray luminosities 

10 8 < Lopt/x < 10 9 L . 

Figure [5] plots the observed optical and X-ray lumi- 
nosities for each SN. For SNe 2008gm and 2008en, the 
optical and/or X-ray luminosities are fairly consistent 
with the contours in Figure |U The optical and/or X- 
ray luminosities may seem a little low compared to the 
expected values given by the dust parameters in Fig- 
ure QJ but we stress that the observed dust parameters 
arc only upper limits due to the limitations of a single- 
component fit. Furthermore, even in an optically thin 
scenario (i.e., r < 1), 50% of the optical flux can be ab- 
sorbed by the dust. Overall, the optical/X-ray fluxes are 
sufficient to power the mid-IR emission. Although Fig- 
ure [2] excl udes previou s ly pub li shed optical data , obser - 
vati ons bvlSmith et al l (|2009bO . [Stritzinger et al.l (|2012t) . 
and lZhang et al.l (J2012 ) yield similar conclusions for SNe 
2005ip, 2006jd, and 2010J1, respectively. Due to the high 
background flux from their underlying galaxies, upper 
limits for SNc 2005cp, 2007rt, and 2008cg (see Table E]) 
do not rule out the presence of high-energy emission from 
CSM interaction. No observations were acquired of SN 
2006qq. 

SN 2010J1 does stand out, however, with a relatively 
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Figure 5. Plots of the Ha line for the spectra in Figure [3] These 
plots show both low- and moderate-resolution spectra obtained 
with Keck LPJS and DEIMOS, respectively. A post-shocked CSM 
will typically have velocities > 100 km s _1 . These relatively broad 
lines are prominent in the LRIS spectra of SNe 2005ip and 2006jd, 
but more ambiguous in the other LRIS spectra. The moderate- 
resolution DEIMOS spectra differentiate between the broader com- 
ponent of CSM interaction (in the case of SN 2008cn) and the nar- 
rower (< 100 km s _1 ) component likely associated with an under- 
lying H II region (SNe 2005cp and 2008gm). The blueshift (relative 
to the rest wavel ength) of the broader component is likely due to 
dust absorption (Smith ct al. 2009b; Stritzingcr ct al. 2012). 



large dust temperatur e for the derived blackbody radius. 
lAndrews et al.l ([201 1[ ) derive similar parameters from 
Spitzer observations obtained on day 90 post-explosion 
and suggest that these characteristics are more consis- 
tent with an IR echo powered by the initial SN flash. 
Optical and X-ray observations, however, reveal clear 
signatures of CSM interaction. Continued monitoring 
will ultimately constrain the dust-shell size and heating 
source. 

Although Spitzer did not detect SNe 2005gn and 
2008J, the dust shell upper limits (0.6 and 0.036 ly, re- 
spectively) derived in §3.11 offer valuable constraints on 
their possible heating mechanism. Assuming dust tem- 
peratures (737 K and 704 K) derived fr o m the first epoch 
of Spitzer observations by iFox et all (|2011l ). Figure [4] 
shows that the dust never required an optical and/or 
X-ray luminosity greater than 10 85 L Q . While we do 
not have late-time optical observations of these two SNc, 



such limits are consistent with the other optical and X- 
ray observations. 

The observed optical and X-ray luminosities are con- 
sistent with the physical scenario described above, but 
the source of the emission is not immediately clear. For 
example, the emission may originate from ongoing CSM 
interaction, but it may also be associated with an under- 
lying H II region. To identify distinct signatures of CSM 
interaction, we obtained ground-based optical spectra of 
the four SNe detected at visual wavelengths (SNe 2005ip, 
2006jd, 2008en, and 2008gm) and SN 2005cp, as summa- 
rized in Table II and Figure El While SN 2005ip exhibits 
clear signs of CSM interaction via the high-ionization Fc 
lines and broadened Ha, the evidence in the other SNe 
is more ambiguous. 

Unlike the LRIS spectra, the moderate-resolution 
DEIMOS spectra resolve down to < 100 km s" 1 and 
can differentiate between H II regions and CSM inter- 
action, which typically has post-shock velocities >100 
km s _1 . Figure [5] illustrates these differences with 
the Ha line. SNe 2005ip and 2006jd both have dis- 
tinct, relatively broad (> 1000 km s _1 ) features asso- 
ciated with the forward shock and a narrower (~ 100- 
200 km s _1 ) c omponent likely associated with the pre- 
shocked CSM (|Fox et alJl201ll ). The line-profile asym- 
metry (a more prominent blue side) of the broader com- 
ponen t is likely due to dus t absorption (|Smith et al.l 
I2009bt iStritzinger et al.ll2012D . The broad feature in SN 
2008en is less obvious, particularly in the LRIS spectrum. 
The moderate-resolution DEIMOS spectrum, however, 
resolves the two components. The LRIS spectrum of SN 
2008gm also docs not rule out a component of ~ 500 
km s^ 1 . The DEIMOS spectra of SNc 2005cp and 
2008gm, however, show only a narrow (< 100 km s _1 ) 
line suggestive of an H II region. 

4. DISCUSSION 

4.1. Mass-Loss History 

The mid-IR traces the characteristics of the CSM at 
the dust-shell radius. Assuming a dust-to-gas mass ratio 
expected in the H-rich envelope of a massive star, Z<j = 
Md/Mg £3 0.01, the dust-shell mass can be tied to the 
progenitor mass-loss rate, 



M- 



Z d Ar l 
3/Md 
4VM Q 



(3) 



for a progenitor wind speed v w . The relatively narrow 
lines of SNe Iln originate in the slow pre-shocked CSM 
and therefore correspond to the progenitor wind speed. 
For these SNe, the narrow-line profiles have a full width 
at half-maximum intensity (FWHM) in the range ~ 100- 
500 km s~ 1 (|Fox et al.l 120 111 ). Assuming a thin shell, 
Ar/r = 1/10, Table [6] lists the associated mass- loss rate 
for each SN. 

By contrast, the optical and/or X-ray observations 
trace the CSM characteristics at the inner radii. As- 
suming a steady mass loss, the rate can be written as a 
function of the optical/X-ray l uminosity, progenitor wind 
speed, and shock ve locity (e.g.. lChugai fc DanzigerllT994l ; 
ISmithet al.ll2009bD : 
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Also plotted are pre-existing mid-IR , optical, and X-ray 

1.I I2012I ; 



Figure 6. Individual plots of the SNe in Figure [2] having the best coverag 

photometry fllmmler & Poolevl [20071: ISmith et al]l2009bl: IFox et al.ll20TH : I Andrews et~aT|[20Tll : IChandra eraLll2012ai rb1: IZhang'et all 
Stritzinger ct al. 2012). While the mid-IR luminosities are comparable, the observations are sparse and comparisons are difficult. The 
optical and X-ray evolution, however, show significant contrasts that likely originate in the progenitor. Multi-wavelength light curves such 
as these offer useful constraints on theoretical models discussed in the text. 



Table 6 

Mass-Loss Rates 



SN M(IR) M(opt/X) 

(M yr" 1 ) (M Q yr" 1 ) 



2005cp 

2005ip 

2006jd 

2006qq 

2007rt 

2008cg 

2008en 

2008gm 

2010J1 



3.4xl0" 2 
8.3xl0" 2 
1.7X10- 1 
1.2X10" 1 
l.OxlO" 1 
6.1xl0™ 2 
7.8xl0- 2 
2.8xl0" 3 
8.8xl0" 2 



1.7xl0" 2 
2.2xl0~ 3 



l.lXlO -3 

1.2xl0" 4 



2v w 
M= — ^--Ew/Xj 

evi ' 



= 2.1 x 10" 



(5) 



J opt/X 



*m 



120 km s" 



3 x 10 41 erg s" 

r)Go^=r)~W', (6) 



where e < 1 is the efficiency of converting shock kinetic 
energy into visual light. While the conversion efficiency 
varies greatly depending on shock speed and wind den- 
sity, we assume e « 0.5, acknowledging that this value 
may be high. Table [6] also lists the mass-loss rates for 
the SNe at the inner radii based on the optical/X-ray 
luminosities in Table [3] 

For the SNe with both IR and optical observations, 
Table [5] highlights that the dust shells correspond to a 
period of higher mass loss and densities. By contrast, the 
CSM that generates the optical and/or X-ray emission 
corresponds to a period of lower mass-loss rates. The 
progenitor mass loss was variable, not continuous. The 
rates given in Table [6] (e.g., M « 10" 1 to 10~ 3 M yr" 1 ) 
are consistent with previous comparisons of SNe Iln to 



the e pisodic dense winds observed in some massive s tars 
(e.g. JFransson et al.llJMJl ISmith et aIll2009all20Tl . In 
almost every case, the dust-shell radius and progenitor 
wind speeds suggest outbursts that occurred on the order 
of tens to hundreds of years prior to the SN explosion. 

4.2. Theoretical Models and Progenitor Diversity 

While stellar evolution models are still limited in 
their ability to achieve such substantial mass loss, the 
nature of the CSM wind can constr ain the progeni- 
tor m ass and explosion physics (e.g., iBalberg fc Loebl 
1201 If ) . The observations in this paper present only a 
single snapshot of the CSM environment, but a num- 
ber of theoretical models can generate the entire multi- 
wavelength evolution as a function of the progenitor 
mass-loss history. For example, the early (< 100 d) op- 
tical and X-ray evolution can be described as a func- 
tion of the immediate CSM density and explosion en- 
ergy (e.g., iChevalier fc Irwinl 120111: iMoriva et all 120111; 



iChevalier fc Irwinl 120121 ; iMorivafc Tominagal 12012ft 



Furthermore. lOfek et afl l|2012ft and lSvirski etall (|2012l ) 
show that for a sufficiently dense CSM, the X-rays will 
escape only months to years after the SN peaks in the 
visible. The escape of the X-rays (sometimes referred 
to as a circumstellar breakout) will di ffer observationally 
for va rying progenitor mass-loss rates (|Chevalier fc Irwml 
|2011[ ). To fully understand the mass- loss history of 
the progenitor therefore requires well-sampled multi- 
wavelength observations. 

Figure [5] plots the multi-wavelength data available for 
individual SNe in our sample, including data from Figure 
[2] and pre-existing optical and X-ray p hotometry for SNe 
2005ip, 2006jd, 2007rt and 2010il fllmmler fc Poole ' 



20071; Smith et al.ll2009b) IFox et al.lfeo 111 IChandra et al 



2012a|bt iZhang et all feoia iStritzinger et all 12012ft 



While fitting the theoreti cal models (such as th ose by 
IMoriva fc Tominagal [20121 and lSvirski et al.ll2012l ) is be- 
yond the scope of this paper, Figure [6] serves_asa poten- 



tially useful template for theorists (e.g. JPan et al.ll2013ft 
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Furthermore, the figure highlights the diversity among 
SNc Iln at various wavelengths. For example, the mid-IR 
luminosities are similar and serve as useful indicators for 
shock interaction and optical/X-ray emission. The op- 
tical and X-ray evolution differ significantly, suggestive 
of different progenitor mass-loss characteristics and, per- 
haps, different progenitors. In the case of SN 2010J1, the 
absorbed luminosity increases by a factor of 2 more than 
1 yr after the visible maximum, revealin g the first case 
of an X -ray breakout from a dense CSM ([Chandra et al.l 
I2012bt l. By contrast, SN 2006jd shows an X-ray plateau 
but no sign of a shock breakout, while SN 2005ip has 
insufficient X-ray observations. The optical light curve 
of SN 2010J1 is quite luminous (> I0 9 L Q ), but fades 
throughout the extent of the observations, whereas SN 
2006jd peaks at ^500 d and SN 2005ip plateaus through 
more than 1000 d. In the case of SN 2006jd, a large X-ray 
luminosity implies a high density, but a small amount 
of p hotoabsorption in the spectrum implies a low den- 
sity ([Chandra et alll2~012aD . These conflicting measure- 
ments suggest an asymmetric CSM. The slow evolution 
of the unabsorbed emission further suggests a progeni- 
tor wind-density profile that does not follow the stan- 
dard r~ 2 model. In the case of SN 2005ip, the late- 
time optical luminosity is over a magnitude larger than 
the X-ray luminosity. SN 2005ip does, however, exhibit 
very high-ionization coronal lines in the visible spectrum. 
iSmith et al.l (j2009bl ) explain these characteristics with 
a lower density CSM interspersed with denser clumps. 
The shock interaction with the dense clumps generates 
X-rays, but with less intensity than in SN 2006jd. Since 
the clumps have a small filling factor, the X-rays can 
escape and photoionize the dense, pre-existing CSM at 
larger radii. 

5. CONCLUSION 

This paper presents multi- wavelength follow-up obser- 
vations of a sample SN Iln population i dentified by late - 
time mid-IR emission from warm dust ([Fox et alj|2011l ). 
For at least five of the nine detected SNe, the data iden- 
tify the predominant heating source of the dust as radia- 
tive optical and/or X-ray emission continuously gener- 
ated by ongoing CSM interaction. Optical s pectra of SNe 
2005ip, 2006jd, 2008en, 2008gm, and 2010J1 (|Smith et al.l 
l2012f ) confirm the CSM interaction with the presence of 
relatively broad Ha lines. The optical spectra of SNe 
2005cp and 2008gm reveal only narrow lines, suggestive 
of an H II region, but these spectra were obtained > 6 
months following the initial Spitzer detection. The IR 
emission may have faded, or radiative shock emission 
does not provide the dominant heating mechanism. Also, 
Spitzer did not detect SNe 2005gn or 2008J, consistent 
with the forward shock overrunning and destroying the 
dust shell. The nondctections place upper limits on the 
dust-shell size, which are also consistent with the radia- 
tive shock model despite the lack of shorter wavelength 
observations. Further Spitzer observations will constrain 
the remaining dust-shell radii and resolve some of the 
ambiguity surrounding SNe 2005cp and 2008gm. 

The multi-wavelength data trace the mass-loss history 
of the progenitor. All of the SNe in this sample appear 
to have undergone an outburst on the order of tens to 
hundreds of years prior to the SN explosion. More re- 
cently, the mass loss appears to be more sporadic. The 



multi-wavelength light curves vary significantly, partic- 
ularly at optical wavelengths, suggestive of a range of 
Type Iln progenitors. However, Spitzer did not detect 
20 of the SNe Iln that occurred between 2003 and 2008. 
From Equation a shock velocity v s rs 10 3 km s _1 into 
a CSM produced by a steady wind with a mass-loss 
rate of only 10 -4 M Q yr -1 emits a luminosity of only 
~ 10 6 Lq, which would be below the detection thresh- 
old for most of these SNe. Furthermore, lower mass- 
loss rates also correspond to lower densities, which result 
in lower dust masses. These nondctections likely corre- 
spond to progenitors that had mass-loss rates < 10 -4 M 
yr -1 , more in line with LBVs in their S Doradus state 



([Humphreys fc Davidson! 1 19941 ). 

The diversity amongst this sample parallels the larger 
SN Iln population. The relatively narrow lines and 
dense CSM associated with SNe Iln have now been 
identified in an unexpectedly divers e list of subclasses , 
including SLSNe (e.g., S N 2006gv: ISmith et all 120071) . 
Type IIP (e.g., SN 2 011ht:lMauerhan et al\\20m Typ e 
IIL (e.g., SN 1998S: ILluet all 120001: IPozzo et al.ll2004D 
Type Ibc (e.g., SN 20 06ic:lFolev et al.|l2007t ISmith et all 
120081) . "impostors" dKochanek et all I2012D. and even 
Type la (e.g. , SN 2 02ic. SN 2008J:|rIamuv et al.l l200l 
iTaddia et all I2012t iSilverman et al.l I2013J )" Existing 
multi- wavelength data do not come close to sampling the 
wide diversity of SNe Iln. Future, well-sampled, multi- 
wavelength observations over year-long timescales will be 
required to classify these different types based on the na- 
ture of the X-ray escape time and the precise relationship 
between CSM interaction and warm dust emission in the 
mid-IR. 
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